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By tuning the radiative coupling of localized surface plasmons to diffracted orders, we demonstrate 
how stop-gaps in plasmonic crystals of nanorods may be opened and tuned. The stop-gap arises from 
the mutual coupling of surface lattice resonances (SLRs), which are collective resonances associated 
with counter-propagating surface polaritons. We present experimental results for three different 
nanorod arrays, where we show how the dispersion of SLRs can be controlled by modifying the size 
of the rods. Combining experiments with numerical simulations, we show how the properties of the 
stop-gap can be tailored by tuning a single structural factor. We find that the central frequency of 
the stop-gap falls quadratically, the frequency width of the stop-gap rises linearly, and the in-plane 
momentum width of the standing waves rises quadratically, as the width of the nanorods increases. 
These relationships hold for a broad range of nanorod widths, including duty cycles of the array 
between 20% to 80%. We discuss the physics in terms of a coupled oscillator analog, which relates 
the tunability of the stop-gaps to the coupling strength of plasmonic modes. 

PACS numbers: 42.70.Qs, 78.67.Qa, 73.20.Mf, 42.25.Fx 



Photonic bandgaps, i.e., frequency regions within 
which optical propagation in a photonic crystal is for- 
bidden, enable to control the flow of light. Gaps arise 
in the dispersion of photonic modes as a consequence of 
the optical periodicity, and the width of the gap is deter- 
mined by the degree of field concentration in the different 
dielectric region^. A more complex scenario arises in 
plasmonic crystals, where light-matter interactions are 
no longer dominated by the optical periodicity alone. 
With the penetration of the electromagnetic field into the 
metal and the excitation of surface plasmon modes, fas- 
cinating phenomena fluorish in plasmonic crystals. For 
this reason, the coupling of surface modes in periodic 
metallic nanostructures has attracted much interest since 
early investigations^, especially for its connection with 
frequency stop-gap^. Coupled surface modes have been 
observed in metallic gratingJ^, subwavelength hole ar- 
ray^^, nanoslit array^, and particle arrays coupled to 
waveguide modes^. Despite the numerous structures that 
have been investigated, very few studies have discussed 
the physical origin of these gaps. One notable excep- 
tion is the work by Barnes and co-workers^, where the 
origin of the gap in metallic sinusoidal gratings is dis- 
cussed. However, the analysis therein contained is not 
easily extended to more complex plasmonic structures, 
which renders difficult the emergence of a simple, intu- 
itive explanation on the origin of the gap. 

It is the aim of this work to provide an intuitive frame- 



work in which the opening of frequency stop-gaps in plas- 
monic crystals can be understood in terms of the cou- 
pling strength of the surface modes involved. We inves- 
tigate nanorod arrays supporting Surface Lattice Reso- 
nances (SLRs), which are dispersive and spectrally nar- 
row collective resonances arising from the diffractive cou- 
pling of localized surface plasmonJ^^^^. This coupling 
occurs near the condition at which a diffraction order 
changes from radiating to evanescent, i.e. at the Rayleigh 
anomaly. In a recent work, we discussed the coupling of 
bright and dark SLRs in nanorod arrays^^. The associ- 
ated stop-gap, modal symmetries, and the very high qual- 
ity factors of SLRs were therein discussed. In this paper, 
we demonstrate how stop-gaps associated with coupled 
SLRs can be selectively opened by tuning the radiative 
coupling of localized surface plasmons to diffracted or- 
ders. We present experimental results for three different 
arrays with varying nanorod dimensions but equal lat- 
tice constants. By combining experiments with numeri- 
cal simulations, we elucidate the influence of the width 
of the nanorod on the dispersion of SLRs. We also find 
scaling laws for the properties of the gap as a function 
of the nanorod width. As we will show, these scaling 
laws are related to the coupling strength of the surface 
modes involved, and they are valid for a wide range of 
duty cycles. 

We have investigated experimentally three 3x3 mm^ 
arrays of gold nanorods fabricated on a silica sub- 
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FIG. 1. (Color online) Experimental extinction spectra as a function of the wave vector component parallel to the surface of 
the array. Scanning electron microscope images of the nanorod arrays yielding the spectra in Figures (a) and (c) are shown 
in the corresponding insets; the scale bar denotes 500 nm. The black solid and dashed lines indicate the (+1,0) and ( — 1,0) 
Rayleigh anomalies, respectively. The incident light is polarized parallel to the width of the nanorods (y-direction) , which is 
(a) 85 nm, (b) 95nm, and (c) 115 nm. The broad, dispersionless resonance on the high frequency side of each spectrum is the 
dipolar localized surface plasmon resonance along the width of the nanorods. The narrower and dispersive resonances below 
the Rayleigh anomalies are the surface lattice resonances. 



strate using Substrate Conformal Imprint Lithography 
(SCIL)^^. The arrays were embedded in a uniform sur- 
rounding medium by placing a silica superstrate preceded 
by n=1.45 index matching fluid to ensure good optical 
contact. The three arrays have lattice constants ax = 600 
nm and ay = 300 nm, but they comprise nanorods which 
differ in size. The nanorods have an approximately rect- 
angular shape in the plane of the array, and a height 
= 38 ± 2 nm. The rod size (length x width) was tuned 
by varying the exposure dose of the electron beam when 
preparing the master for nanoimprint. This procedure 
yielded rods of size (a)420 x 85 nm^, (b)430 x 95 nm^, 
and (c)450 x 115 nm^, which correspond to the measure- 
ments in Figures [TJa), [ijb), andjljc), respectively. The 
tolerances of these in-plane dimensions are on the order of 
±10 nm. Figures TJa) and[TJc) show a Scanning Electron 
Microscope (SEM) image of the corresponding array, and 
a cartesian triad in the inset which we use to describe the 
measurements next. 

Figure [l] shows the extinction, defined as 1 — T with 
T the zeroth order transmittance, for the three arrays 
described above. The extinction is displayed as a func- 
tion of the reduced frequency, i.e., the angular frequency 
normalized by the speed of light in vacuum, and the 
component of the incident wave vector parallel to the 
surface, which is given by k\\ = ^nsin(^)x, with n the 
refractive index of silica and the angle of incidence. 
The sample was rotated around the y-axis while the y- 
polarized collimated beam from a halogen lamp impinged 
onto the sample, probing the short axis of the nanorods. 
The broad, dispersionless extinction peak seen on the 
high frequency side of the spectra for all three arrays 
corresponds to the excitation of Localized Surface Plas- 
mon Resonances (LSPRs) in the individual nanorods. 
The black solid and dashed lines indicate the (+1,0) 
and (-1,0) Rayleigh anomalies of the arrays, respectively. 



The Rayleigh anomalies are solutions to the equation 
klut = ^In sm2(<^) + ml{2Ti/axY + m^(27r/a^)2 + 2kin 
sin(6>)mi(27r/aa^), where kin and kout are the modulus of 
the incident and scattered wave vectors and rrij (j = 1, 2) 
are the integers defining the order of diffraction. Phys- 
ically, these so-called anomalies represent the frequency 
and wave vector for which the corresponding diffracted 
orders are propagating grazing to the surface of the ar- 
ray. The coupling of LSPRs to the Rayleigh anomalies 
gives rise to the SLRs, which manifest in the measure- 
ments as narrow and dispersive peaks in extinction at 
slightly lower frequencies than the associated Rayleigh 
anomalies. 

Figures [TJa)-(c) show the gradual opening of a stop- 
gap in the dispersion relation of SLRs as the nanorod 
width increases. The gap is centered near 7 mrad/nm 
in Fig. [ija), but its central frequency is lowered and its 
width Augap increases as the nanorods become wider. 
We note that this is not a complete photonic band-gap, 
since it only exists for light polarized parallel to the short 
axis of the nanorods and with an in-plane wave vector 
component parallel to the long axis of the nanorods. For 
light polarized parallel to the long axis of the nanorods, 
the dipolar LSPR lies at lower energies than the (±1, 0) 
diffraction orders at normal incidence, which results in a 
weak diffractive couplin^^. On the other hand, for an in- 
plane wave vector component parallel to the short axis of 
the nanorods, the (±1,0) Rayleigh anomalies are degen- 
erate, leading to degenerate (±1,0) SLRs and therefore 
no gapP^. 

Inspired by previous work explaining electroma gnetic 
resonance phenomena in terms of coupled oscillator^^, 
we recently introduced an analog to the plasmonic crys- 
tal consisting of three mutually coupled harmonic oscil- 
lators^^. In this analogy, the conduction electrons in the 
nanorod driven by the electromagnetic field are modeled 
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FIG. 2. (Color online) The black open circles, blue open dia- 
monds, and red crosses, are the experimental extinction spec- 
tra at normal incidence for arrays of nanorods of width 85 
nm, 95nm, and 115 nm, respectively. The black solid, blue 
dashed, and red dash-dot curves are calculations of the ab- 
sorbed power in the coupled oscillator model described in the 
text and with fitting parameters as given in Figure 3(c). The 
dashed vertical line indicates the eigenfrequency UJ2 = 7.3 
mrad/nm representing the Rayleigh anomaly in the oscillator 
model. 

as oscillator 1 driven by a harmonic force F = Foe"'^'^^^. 
The (+1,0) and (-1,0) Rayleigh anomalies are modeled by 
oscillators 2 and 3, respectively. The equations of motion 
of the coupled system are, 

Xi + 7iXi + Uj\xi - ^12^2 - ^13^3 = F, 

X2 + 72^2 + ^2^2 - ^12^1 - ^23^3 = 0, (1) 

Xs + 73^3 + ^3^3 - ^13^1 - ^23^2 = 0, 

with Xj, 7j , and ujj {j = 1,2,3) the displacement from 
equilibrium position, damping, and eigenfrequency as- 
sociated with the j^^ oscillator, respectively, and Qjk 
{k = 1,2,3 and j ^ k) the coupling frequency between 
the j^^ and k^^ oscillator. The absorbed mechanical 
power by oscillator 1 from the driving force, given by 
P{t) = Fxi^ is representative of the extinct optical power 
by the array. Integrating P{t) over one period of oscil- 
lation and scanning the driving frequency Ug yields an 
absorbed power spectrum, which we compare to the ex- 
tinction measurements. 

Figure |2] displays the extinction spectra for the three 
arrays at normal incidence, together with calculations 
done with the coupled oscillator model as fits to the 
measurements. Only the (+1,0) SLR appears in the 
spectra, so the model reduces to two oscillators, i.e., 
^13 = ^23 = 0. The third oscillator is uncoupled be- 
cause the (-1,0) SLR is a dark state at normal incidence. 
This condition arises due to an antisymmetric (quadrupo- 
lar) character of the (-1,0) modal fields, which results in 



a vanishing extinction at normal incidence^^. Figure [ij 
shows that for an increasing nanorod width the LSPR 
displays a diminishing center frequency and a broader 
linewidth. This behavior is associated with increased re- 
tardation and radiative damping^^, which we model in 
the calculations of Fig. [2] by lowering ui and increasing 
7i. The eigenfrequency of oscillator 2, representing the 
Rayleigh anomaly frequency, is set to co'2 = 7.3 mrad/nm 
for all three cases. In order to unambiguously demon- 
strate the critical role of the coupling between LSPRs 
and the Rayleigh anomaly in determining the SLR line- 
shape, we set the damping of oscillator 2 (whose coupling 
gives rise to the SLR) equal in all three cases, given by 
72 = 0.001 mrad/nm. Considering that the losses are ex- 
pected to increase for the bigger nanorods, this is unlikely 
the exact case in the experiments. However, as we show 
next, fixing 72 allows us to see the effect of changing the 
coupling frequency 1^12, which stands for the radiative 
coupling strength between the LSPR and the Rayleigh 
anomaly. 

The measurements in Fig. [2] show two main effects on 
the SLR lineshape as the width of the nanorods increases: 
i) the peak resonance frequency is increasingly detuned 
from the Rayleigh anomaly, and ii) the linewidth broad- 
ens. As it is shown next, this behavior is explained by the 
change of the coupling constant between the LSPR and 
the Rayleigh anomaly. In Fig. [3]^ a) we plot the detuning 
of the SLR from the Rayleigh anomaly, oora—^slr^ and 
the linewidth at Full- Width Half Maximum (FWHM) of 
the SLR as a function of the width of the nanorods. It is 
remarkable that these two quantities are in quantitative 
agreement, which indicates a direct connection between 
radiative losses and the peak resonance frequency with 
respect to the corresponding Rayleigh anomaly. The im- 
plications of this connection for sensing small changes to 
the bulk refractive index by means of plasmonic nanopar- 
ticle arrays have been recently discussed^. A universal 
scaling of the figure of merit of plasmonic sensors, which 
is a function of the detuning ujra — ^slr alone, was 
therein found. In Fig. |3jb) we plot ujra — ojlspr as 
a function of the nanorod width, which displays the di- 
minishing frequency difference between the LSPR and 
the Rayleigh anomaly. This decrease in the magnitude 
of ura — ojlspr with a simultaneous broadening of the 
LSPR promotes a stronger radiative coupling of localized 
surface plasmons to diffracted orders, i.e., an increase in 
Figure [3|c) shows the values of 1^12 used to fit the 
measurements in Fig. [2] Although in the fitted range 1^12 
may well be described by a linear function, we have used 
a quadratic function for a reason that will be clarified fur- 
ther in the text. We see that increasing 1^12 detunes the 
SLR from the Rayleigh anomaly and broadens it in the 
right amount to have an excellent agreement with the 
measurements. A crucial understanding in attributing 
the observed behavior to 1^12 mainly lies in the fact that 
changing 72 (the intrinsic damping of oscillator 2) can not 
lead to the right detuning ujra — ^slr- We have veri- 
fied this through several calculations (not shown here) 
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FIG. 3. (Color online) Figure (a) shows the (+1,0) SLR 
linewidth at Full- Width Half Maximum (FWHM) as blue 
open circles, and the detuning of the SLR from the Rayleigh 
anomaly, ura — ojslr^ as red open triangles. Both quanti- 
ties are taken from the extinction measurements at normal 
incidence. The solid and dashed lines are guides to eye for 
the FWHM and detuning, respectively. Figure (b) shows the 
frequency difference between the Rayleigh anomaly and the 
Localized Surface Plasmon Resonance (LSPR), ujra—ojlspr^ 
as green open diamonds. The green dashed line is a guide to 
eye. Figure (c) shows the coupling frequency Q12 used in 
the coupled oscillator model to calculate the spectra shown 
in Figure 2. The black dashed curve is a quadratic fit to Q12. 
All quantities in (a) and (b) are plotted as a function of the 
nanorod width. The error bars in (a) and (b) stem from the 
resolution in the spectrometer, whereas the error bars in (c) 
represent the uncertainty in the fitting. 



in which 72 and Vti2 were varied independently and/or 
simultaneously. Next, we consider the more general case 
of inclined incidence light, where the three oscillators are 
mutually coupled. 

Inspection of Fig. [l] points towards the connection be- 
tween the formation of standing waves in the high fre- 
quency SLR band and the opening of the gap. We con- 
sider the gap to open at the slowdown point for the high 
frequency band, i.e., where the dispersion of the (+1,0) 
SLR flattens and the group velocity slows down to zero. 
Figure [4] shows the extinction of the three arrays at the 
slowdown point, which occurs at (a) k\\ = 0.13 mrad/nm, 
(b) k\\ = 0.18 mrad/nm, and (c) k\\ = 0.35 mrad/nm. 
In the same Figure we plot the absorbed power spec- 
tra for the oscillator model fitting each measurement. 
From the frequency and linewidth of the dispersionless 
LSPR we set the eigenfrequency and damping of oscilla- 
tor 1 to (a) uji = 9.5 mrad/nm, 71 = 2.3 mrad/nm, (b) 
uji = 9.2 mrad/nm, 71 = 2.4 mrad/nm, and (c) uji = 8.8 
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FIG. 4. (Color online) The black open circles are cuts of the 
measured extinction spectra for the three arrays in Figure 1, 
following the same labeling convention for the antenna width, 
i.e., (a) 85 nm, (b) 95 nm, and (c) 115 nm. For each case 
the extinction is shown at the slowdown value of , which 
is at (a) k\\ = 0.13 mrad/nm, (b) = 0.18 mrad/nm, and 
(c) k\\ = 0.35 mrad/nm. The green solid curves represent the 
absorbed power in the coupled oscillator model described in 
the text, with coupling and damping frequencies as given in 
Table 1. The blue and red dashed lines are the (+1,0) and 
(-1,0) Rayleigh anomalies as predicted by the theory (conser- 
vation of the parallel component of the wave vector across the 
interface), respectively. The blue and red dash-dot lines are 
the eigenfrequencies 002 and 002, , respectively, set in the oscilla- 
tor model to calculate the absorbed power spectrum for each 
case. 



mrad/nm, 71 = 2.7 mrad/nm. The coupling and damp- 
ing frequencies used to fit the SLR lineshapes are given 
in Table 1, together with the eigenfrequencies 002 and uj^ 
representing the Rayleigh anomalies. The blue and red 
dashed vertical lines in Fig. [4] indicate the (+1,0) and (- 
1,0) Rayleigh anomalies, respectively, as predicted by the 
equation describing the conservation of the wave vector 
component parallel to the surface of the grating. These 
are the frequencies at which the black lines in the disper- 
sion diagrams in Fig. [l] cross the value of k\\ inspected 
in cases (a)-(c). The blue and red dash-dot vertical lines 
in Fig. |4] indicate where the dips in extinction associated 
with the Rayleigh anomalies are seen in the experiment, 
which are also the eigenfrequencies used in the oscillator 
model. We see in Fig. |4] that as the nanorod widens and 
the gap opens there is an increased frequency deviation of 
the extinction dip with respect to the corresponding theo- 
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FIG. 5. (Color online) Figure (a) is a zoom into the stop- 
gap displayed in the measurements of Figure 1(a). Figure (b) 
shows FDTD simulations results for an array of comparable 
antennas. The extinction is displayed by the same color scale 
for both measurements and simulations. The black solid and 
dashed lines represent the (+1,0) and (-1,0) Rayleigh anoma- 
lies, respectively. 



retically predicted Rayleigh anomalies. Figure |2] displays 
the same phenomenon for normal incidence light, but 
the effect is much less pronounced than at the slowdown 
point. A small shift between the diffraction edge and 
the associated extinction dip was also observed by Augie 
and Barnes in the normal incidence extinction spectra of 
various metallic nanoparticle arrays^^, but its origin was 
not discussed. It is crucial to realize that when deriving 
the conditions for which a diffracted order propagates 
grazing to the surface of the grating, i.e. the Rayleigh 
anomaly equation yielding the black lines in Fig. [l] and 
the dashed lines in Figs. [2] and [4j the form factor of the 
grating is not considered. The derivation follows from 
equating the incoming and outgoing wave vectors, the 
former including the momentum added by the grating. 
However, this first order analysis neglects any changes in 
momentum that could arise from the dimensions of the 
particles constituting the grating. Since the dimensions 
of the individual particles determine their polarizability, 
i.e., the frequency and damping of the associated LSPR, 
we propose that the small shift in frequency of the extinc- 
tion dip with respect to the corresponding theoretically 
predicted Rayleigh anomaly may be rooted in the cou- 
pling strength of LSPRs to diffracted orders. Indeed, we 
observe in Figs, [l] |2] and [4] that for the widest nanorods, 
where the coupling strength of the LSPR to the (± 1,0) 
Rayleigh anomalies is highest, the aforementioned shift 
is also highest. 

In order to investigate the properties of the gap for 
a broader range of nanorod widths, we have performed 
Finite Difference in Time Domain (FDTD) simulations 
with an in- house developed software. We have calculated 
dispersion relations in extinction for arrays with lattice 
constants ax = 600 nm and ay = 300 nm, and rods of 
dimensions L = 450 nm, h = 40 nm. The width of the 
rods was varied between w = 20 nm and w = 280 nm, in 
steps of 20 nm. Although the three experimental arrays 
displayed a small change in the length of the nanorods, 
the dominant contribution to the properties of the gap 



stems from the width of the nanorod (within the range of 
length variation). This assumption is based on the fact 
that the coupled surface modes herein discussed arise for 
light polarized parallel to the width of the nanorods, and 
the associated LSPR red-shifts as the width increases, 
which is the expected behavior due to the depolarization 
field along this dimensioil^. We validate the simulation 
results in Fig. |5j where we compare the measured and 
calculated dispersion relations near the gap for an array 
of gold nanorods of 80 nm width in both cases. From the 
dimensions of the three experimental arrays previously 
given, it can be recognized that the nanorods in the sim- 
ulations are about 7% shorter in length. Nevertheless, 
a good qualitative agreement is observed between mea- 
surements and simulations. 

In Fig. [6]^a) we plot the frequency width of the stop- 
gap, Augap^ and the k\\ width of the standing waves in 
the high frequency SLR band, Aksw^ both as function 
of the width of the nanorods. Aujgap is taken to be the 
frequency difference between the slowdown point and the 
(-1,0) SLR, both evaluated at the k\\ value of the slow- 
down point. Aksw is equal to the value of k\\ at the 
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FIG. 6. (Color online) Figure (a) shows the frequency width 
of the stop-gap, Augap, and the k\\ width of the stand- 
ing waves in the high frequency SLR band, Aksw- The 
black open squares and solid line represent measurements and 
FDTD simulations results, respectively, for Aujgap. The green 
open circles and dashed curve represent measurements and 
FDTD simulations results, respectively, for Aksw- Figure (b) 
shows the detuning of the (±1,0) SLRs from their respective 
Rayleigh anomalies. The (-1,0) detuning is plotted with red 
left-pointing triangles and a solid line, for measurements and 
FDTD simulations results, respectively. The (+1,0) detun- 
ing is plotted with blue right-pointing triangles and a dashed 
curve, for measurements and FDTD simulations results, re- 
spectively. All quantities in (a) and (b) are plotted as a func- 
tion of the nanorod width. 
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TABLE I. Coupling and damping frequencies used in equation 1 to reproduce the spectra shown in Figure 4. In the entries for 
which a minimum estimate is given, the value in parenthesis is the value used in the model yielding the spectra in Figure 8. 
All quantities are given in units of mrad/nm 
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slowdown point. In Fig.[6[b) we plot the detuning of the 
(±1, 0) SLRs from their respective theoretically predicted 
Rayleigh anomalies, also as a function of the nanorod 
width. The open data points are taken from the mea- 
surements of the three arrays, and the solid and dashed 
curves result from the simulations. A central finding is 
that both Augap and the (-1,0) detuning are linearly in- 
creasing functions of the nanorod width, while Aksw 
and the (±1,0) detuning are both quadratically increas- 
ing functions of the nanorod width. We therefore find 
the frequency width of the gap to be correlated with the 
(-1,0) detuning, and the k\\ width of the standing waves 
to be correlated with the (±1,0) detuning. In terms of 
the coupled oscillator model, this can be interpreted as 
follows: the gap opens with a linear increase in (cou- 
pling of LSPR to the (-1,0) order), whereas the in-plane 
momentum width of the standing waves broadens with 
a quadratic increase in 1^12 (coupling of LSPR to the 
(±1,0) order). The latter dependance is the reason for 
which Qi2 was fitted with a quadratic function in the nor- 
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FIG. 7. (Color online) Figure (a) shows the central frequency 
of the stop-gap, uogap, with black open squares representing 
measurements and a black solid curve representing FDTD 
simulation results. Figure (b) shows the gap-midgap ratio, 
with blue open circles representing measurements and a blue 
solid curve representing FDTD simulations results. All quan- 
tities in (a) and (b) are plotted as a function of the nanorod 
width. 



mal incidence spectra, i.e., the dashed curve in Fig. [3[c). 

Further quantifying the properties of the gap, in 
Fig. [7|^a) we plot the central frequency of the stop-gap, 
^gap^ which falls quadratically for increasing nanorod 
width. In Fig. ^h) we plot the gap-midgap ra- 
tio, Aujgap/ujgap which rises quadratically for increas- 
ing nanorod width. Figure [7|^a) shows that for nanorods 
wider than those considered in the experiments, the cen- 
tral frequency of the gap falls into the infrared part of the 
spectrum. Limited by the spectral response of our acqui- 
sition system (based on a silicon detector), we are unable 
to perform measurements for the wider nanorods. Nev- 
ertheless, in light of the good agreement between mea- 
surements and simulations we believe Figs. [6]^a) and [t] 
accurately encompass all properties of the gap as function 
of the nanorod width, and provide a suitable recipe upon 
which stop gaps in plasmonic crystals may be opened 
and tuned. Comparing to previous work, where it was 
found that the gap's width is a linearly increasing func- 
tion of the modulation amplitude in metallic sinusoidal 
gratings^, we have also found a linear dependance of the 
gap's width on a single structural factor. However, the 
connection between these two works is not an obvious 
one, since whereas Barnes and co-workers have varied a 
dimension out of the plane of propagation, we have inves- 
tigated structures of equal height but variable dimensions 
in the plane of propagation. From a fabrication point of 
view, a precise in plane structuring of plasmonic crystals 
may offer a higher degree of versatility in how stop-gaps 
may be tuned and at a greater ease. 

We now discuss the validity range of the previously dis- 
cussed SLR properties on the nanorod width. Figure [S] 
shows the extinction of two arrays with lattice constants 
ax = 600 nm, ay = 300 nm, nanorods of length= 450 
nm and height = 40 nm, and embedded in a homoge- 
neous environment of n=1.45; these conditions are the 
same as in the experiments presented in Fig. [l] The 
width of the nanorods in Fig. [s] is (a) 60 nm, and (b) 
240 nm, which correspond to the extremes of the curves 
plotted in Figs. [6] and [T] Due to the very different ex- 
tinction efficiency of the two arrays in the spectral region 
of interest, the extinction is displayed by a logarithmic 
color scale common to both graphs, which allows quan- 
titative comparison of the spectra. The low extinction 
and remarkably narrow resonances seen in Fig. [sj^a) can 
be explained in terms of the coupled oscillator model as 
a consequence of the low coupling strength of LSPRs to 
diffracted orders. This point can be inferred from the 
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FIG. 8. (Color online) FDTD simulations results for the 
extinction of an array of gold nanorods of width (a) 60 nm 
and (b) 240 nm. Both arrays have nanorods of length = 450 
nm, height = 40 nm, in a lattice with constants ax — 600 nm 
and tty = 300 nm. The black solid and dashed lines indicate 
the (+1,0) and (-1,0) Rayleigh anomalies, respectively. The 
extinction is displayed by a logarithmic color scale common 
to both graphs. 

quantities given in Table 1 , where all coupling frequencies 
ftjk are seen to decrease as the nanorod width decreases. 
Lower fti2 and ftis therefore translate into a lower cou- 
pling strength of LSPRs to diffracted orders, whereas a 
lower Q23 translates into a lower coupling between the 
(+1,0) and (-1,0) SLRs and therefore a smaller gap. As 
the width of the nanorods increases and the Qjk terms in- 
crease, the extinction at the SLRs also increases, the reso- 
nance linewidth broadens, and the frequency gap widens, 
eventually reaching the case displayed in Fig. [sjb) for a 
nanorod width of 240 nm. Figure [sj^b) displays a dipo- 
lar LSPR near 6.8 mrad/nm, which is lower in frequency 
than the diffraction edge at normal incidence. The (+1,0) 
SLR can still be recognized in the spectrum from the 
non-dispersive feature near 6.0 mrad/nm (in the red tail 
of the LSPR). However, for wider nanorods the LSPR 
shifts to lower frequencies, thereby making the (+1,0) 
SLR indistinguishable from the LSPR. As the energy of 



the dipolar LSPR becomes substantially lower that the 
(±1,0) diffraction orders, diffractive coupling of dipolar 
LSPRs becomes very weak, and the properties of coupled 
SLRs leading to the opening of the stop-gap significantly 
deviate. 

In conclusion, we have shown how the radiative cou- 
pling strength of localized surface plasmons to diffracted 
orders in periodic arrays of nanorods can be tuned. This 
tuning is achieved experimentally by modifying the width 
of the nanorods, and results in the opening of frequency 
stop-gaps whose properties are therefore tunable. A 
quadratic dependence of both the frequency width of the 
stop-gap and the (-1,0) SLR detuning was found on the 
nanorod width. A linear dependence of both the in-plane 
momentum width of the standing waves in the high- 
frequency SLR band and the (+1,0) SLR detuning was 
found on the nanorod width. In light of a coupled oscil- 
lator analog to the plasmonic crystal, we have associated 
these two correlations with the coupling strength of local- 
ized surface plasmons to the (-1,0) and (+1,0) diffracted 
orders. Supporting experiments with numerical simula- 
tions for a wider set of nanorod widths, we have analyzed 
the properties of the gap and discussed the limiting cases 
where the scaling laws we present deviate. Although we 
have only considered nanorod arrays in this work, similar 
results are expected to hold in periodic arrays of particles 
with different geometries whereby diffractive coupling of 
localized surface plasmons is possible. Our results there- 
fore pave the road towards nanoscale light manipulation 
in 2D plasmonic crystals, since stop-gaps allow to selec- 
tively enhance or suppress light-matter interactions. 
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